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ABSTRACT 

Using an atmospheric general circulation model (the National Center for Atmospheric Research Community 
Climate Model: CCM2), the effects on climate sensitivity of several different cloud radiation parametenzations 
have been investigated. In addition to the original cloud radiation scheme of CCM2, four parametenzations 
incorporating prognostic cloud water were tested: one version with prescribed cloud radiative properties and 
three other versions with interactive cloud radiative properties. The authors’ numerical expenments employ 
perpetual July integrations driven by globally constant sea surface temperature forcings of two degrees, both 

positive and negative. . . M 

A diagnostic radiation calculation has been applied to investigate the partial contributions of high, middle 
and low cloud to the total cloud radiative forcing, as well as the contributions of water vapor temperature, and 
cloud to the net climate feedback. The high cloud net radiative forcing is positive, and the middle and low cloud 
net radiative forcings are negative. The total net cloud forcing is negative in all of the model versions. The 
effect of interactive cloud radiative properties on global climate sensitivity is significant. The net cloud radiative 
feedbacks consist of quite different shortwave and longwave components between the schemes with interactive 
cloud radiative properties and the schemes with specified properties. The increase in cloud water content in the 
warmer climate leads to optically thicker middle- and low-level clouds and m turn to negative shortwave feedbacks 
for the interactive radiative schemes, while the decrease in cloud amount simply produces a positive shortwave 
feedback for the schemes with a specified cloud water path. For the longwave feedbacks, the decrease in high 
effective cloudiness for the schemes without interactive radiative properties leads to a negative feedback, while 
for the other cases, the longwave feedback is positive. . . 

These cloud radiation parametenzations are empirically validated by using a single-column diagnostic model, 
together with measurements from the Atmospheric Radiation Measurement program and from the Tropical Ocean 
Global Atmosphere Combined Ocean-Atmosphere Response Experiment. The inclusion of prognostic cloud 
water produces a notable improvement in the realism of the parameterization, as judged by these observations. 
Furthermore, the observational evidence suggests that deriving cloud radiative properties from cloud water 
content and microphysical characteristics is a promising route to further improvement. 


1. Introduction 

The net effect of clouds on the radiation balance of 
the earth, referred to as the cloud radiative forcing, is 
negative (Ramanathan et al. 1989) and has an average 
magnitude of about 10-20 W nr 2 . It consists of a short- 
wave cooling (the albedo effect) of about 40-50 W nr 2 
and a longwave warming (the greenhouse effect) of 
about 30 W nr 2 . The size of the observed average net 
cloud forcing is several times the expected value of the 
direct radiative forcing from a doubling of C0 2 (about 
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4 W m“ 2 ), suggesting that cloud effects can play a sig- 
nificant role in climate change. 

Parameterization of the effects of cloud remains one 
of the main areas of uncertainty in the numerical sim- 
ulation of climate change. The computed changes in 
cloud amount and consequent cloud radiative feedbacks 
for a warmer climate have been investigated by many 
authors using general circulation models (GCMs) and 
other tools. In general, GCM results show that cloud 
processes play a key role in determining the sensitivity 
of climate to an external forcing. In comparing 19 
GCMs, Cess et al. (1990) found that the global-average 
climate sensitivity of the models varied by about a factor 
of 3, when measured by surface temperature changes, 
and that much of this difference was related to differ- 
ences in cloud radiation parameterizations among the 
models. Senior and Mitchell (1993) tested four different 
cloud schemes with a version of the U.K. Meteorolog- 
ical Office atmospheric GCM. Schemes that included 
prognostic cloud water variables were more realistic in 
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some respects than schemes that did not, although all 
model versions produced control climates that still 
showed marked differences from the observational data. 
The sensitivity of each of the versions to a doubling of 
C0 2 was critically dependent on the scheme; when cloud 
was parameterized on relative humidity, a strong pos- 
itive cloud radiation feedback was found, while when 
interactive cloud water was included, a negative feed- 
back was found. 

Furthermore, Senior and Mitchell (1993) found that 
there were substantial differences between the results of 
C0 2 doubling experiments and those of inverse climate 
change experiments forced by prescribed changes in sea 
surface temperature. One motivation for the present 
study is to examine the model dependence of these re- 
sults. If we were to test the same cloud schemes in a 
different GCM, the results might be different, because 
of differences in physical parameterizations other than 
cloud processes, for example. Additionally, nonlinear 
interactions between such processes might cause sig- 
nificant differences in the cloud distributions and the 
resulting radiation fields. 

In the present study, we test five different cloud ra- 
diation schemes in the National Center for Atmospheric 
Research (NCAR) Community Climate Model Version 
2 (CCM2). The schemes consist of the original CCM2 
scheme (CCM2), a prognostic cloud water scheme 
(CW), a scheme with explicit cloud water path and fixed 
cloud droplet effective radius (CWRF), one with vari- 
able droplet effective radius for warm cloud only 
(CWRV), and one with variable droplet effective radius 
including ice phase cloud (CWRI). Brief descriptions 
of each scheme are presented in section 2. In section 3, 
a diagnostic method proposed by Zhang et al. (1994) 
using an off-line radiation calculation is briefly de- 
scribed. With this approach, we can examine the role 
of various physical processes in determining the total 
climate feedback. To investigate the sensitivity of each 
cloud radiation scheme, a ±2 K globally constant sea 
surface temperature (SST) perturbation, in conjunction 
with a perpetual July simulation, is adopted as a sur- 
rogate climate change (Cess and Potter 1988; Cess et 
al. 1990; Zhang et al. 1994). As Cess and colleagues 
have pointed out, this procedure is in essence an inverse 
climate change simulation: rather than introducing an 
external forcing into the model — for example, doubling 
C0 2 — and then letting the climate respond to this forc- 
ing, we instead prescribe the climate change and let the 
models in turn produce their respective forcings. To fo- 
cus solely on atmospheric feedback processes, we have 
excluded the ice-albedo feedback by fixing the ice- 
lines — that is, SST in regions adjacent to sea ice is fixed 
at values appropriate to the present July condition — 
disregarding the ±2 K SST perturbation, as in the stud- 
ies cited above. 

Section 4 describes the model simulations of cloud 
and radiation fields of both the present climate and the 
SST-induced climate change, together with satellite- 


derived observational estimates of cloud radiative forc- 
ing fields from the Earth Radiation Budget Experiment 
(ERBE, e.g., Barkstrom 1984) narrow-field-of-view 
(NFOV) monthly mean radiation budget data (this is the 
so-called ERBE S-4 data). We set the solar declination 
and SST at mid-July conditions and have generated a 
240-day perpetual July integration for each case. The 
final 60 days of each integration are examined for this 
study. We have fixed the temperature in the bottom two 
layers of the four-layer finite-heat-capacity soil model 
of the CCM2 at the model’s July climatological values, 
regardless of the SST setting, because the relaxation 
time of the bottom layer is about 5 yr (Zhang et al. 
1994). This allows quick equilibration between the sur- 
face and the atmosphere. However, we retained the di- 
urnal variation of solar insolation. The CCM2-simulated 
present climate shows good qualitative agreement with 
the observational data, but quantitatively there naturally 
are significant disagreements. 

In section 5, we compare the cloud radiation param- 
eterizations using a single-column diagnostic model, to- 
gether with measurements from the Atmospheric Ra- 
diation Measurement (ARM) program and from the 
Tropical Ocean Global Atmosphere Combined Ocean- 
Atmosphere Response Experiment (TOGA COARE). 
Our conclusions are presented in section 6. 

An overall decrease in cloud amount in a simulated 
warmer climate, which is typical in GCM results, can 
significantly reduce the infrared trapping (greenhouse 
effect). Hence, the earth can emit longwave radiation 
more efficiently; that is, the cloud change moderates the 
global wanning (a negative feedback). But, at the same 
time, the incoming shortwave radiative forcings of the 
surface-atmosphere system may increase, because the 
decreases in cloud amount cause a reduction in solar 
radiation reflected by clouds (a positive feedback). Fur- 
thermore, cloud height feedbacks also produce a posi- 
tive feedback in many models. If. in addition, we con- 
sider changes in the radiative properties of clouds, the 
above feedback mechanisms may become much more 
complicated. In this study, we have tested the effect of 
interactive radiative properties, especially microphysi- 
cal cloud radiative properties, on the simulated climate 
change. Somerville and Renter (1984) hypothesized that 
a warmer atmosphere would lead to increased average 
cloud water content, even if there is little change in 
average cloud cover. Using a simple radiative-convec- 
tive model, they found that a warming produced opti- 
cally thicker and so more reflective cloud, tending to 
cool the surface-atmosphere system by providing a neg- 
ative shortwave feedback on the climate change. In this 
study we have tested this hypothesis using a GCM and 
have obtained results consistent with those of Somer- 
ville and Remer. 

2. Model description 

a. Original CCM2 scheme (CCM2) 

There are two standard versions of CCM2 in terms 
of horizontal resolution: a T42 version (approximately 
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equivalent to 2.8° lat X 2.8° long) and an R15 version 
(4.5° lat X 7.5° long). In the vertical, there are 18 levels 
in both the T42 and the R15 versions, with a top at 
2.917 mb. In this study, we use the lower resolution 
version (R 15) for economy. 

The scheme for predicting cloud fraction in CCM2 
is based on the algorithm used in the operational 
ECMWF model (Slingo 1987), to which some modifi- 
cations have been made (Hack et al. 1993). Layer cloud 
fraction depends on relative humidity, vertical motion 
and static stability. Some of the major changes from the 
Slingo (1987) formulation are the following: clouds are 
allowed to form in any tropospheric layer, except the 
layer nearest the surface; low-level clouds do not form 
unless the vertical velocity is upward; and the relative 
humidity thresholds for mid- and upper-level clouds are 
functions of atmospheric stability. However, we have 
eliminated the constraint on cloud-top height of the orig- 
inal CCM2 scheme. 

The convective cloud water content is explicitly es- 
timated in the CCM2 algorithm, during the course of 
computing convective rain rate, but there is no explicit 
cloud water estimation for layer clouds (Hack 1994). 
Therefore, the radiative properties of clouds, for both 
solar and terrestrial radiation, are calculated separately 
in CCM2 using cloud water paths (CWP), which are 
prescribed as a function of latitude, cloud height, and 
cloud-layer thickness. 

The computation of large-scale precipitation in 
CCM2 is based on an adjustment process. If the at- 
mosphere is stable but supersaturated, the moisture field 
is adjusted to be saturated, and the temperature field is 
simultaneously adjusted to reflect the heating due to the 
release of latent heat. The corresponding stable precip- 
itation rate is computed accordingly. See Hack et al. 
(1993) for details. 

b. Cloud water scheme ( CW) 

The cloud water scheme is based on that of Smith 
(1990), as is that of Senior and Mitchell (1993). The 
scheme introduces an additional prognostic variable rep- 
resenting layer cloud water content (q c ). Then, 

9 h - = q + 4c* 

where and q are the total water content and specific 
humidity, respectively. Here, q w is increased by con- 
vergence of water vapor due to convection or large- 
scale advection and by evaporation of precipitation, and 
it is reduced by precipitation. In this study, the large- 
scale advection and the boundary layer vertical mixing 
of q t . are ignored for simplicity. 

The calculation of cloud amount and cloud water con- 
tent uses cloud distribution concepts; that is. a statistical 
distribution of cloud-conserved variables about their 
grid-box mean is assumed. In the mean, the amount of 
cloud and its water content clearly depend on the dif- 
ference between q„ and the saturation specific humidity, 


while the statistics of the local deviation from the mean 
or sub-grid-scale fluctuation, say 5, are defined by a 
prescribed distribution function, say G(s). The mean 
cloud fraction is the fraction of the grid box with cloud 
water content greater than zero, and the mean cloud 
water is the integral of the local cloud water content. 
For this study we adopt G(s) as a simple symmetric 
triangular function, as used by Smith (1990). Unlike 
models without a prognostic cloud water variable, in 
which all condensed water is precipitated, a model with 
cloud water requires a parameterization of the depletion 
of cloud water due to precipitation. We use a scheme 
based on that of Sundqvist (1981). but frozen water is 
assumed to precipitate with a specified fall speed as soon 
as it forms. 

The source code of this scheme is extracted from the 
U.K. Meteorological Office Unified Model. All the pa- 
rameter values are taken to be identical to those in the 
Unified Model, except the values of critical relative hu- 
midity for determining cloud fraction and the cloud wa- 
ter content. These critical values were retuned for re- 
alism of the control climate from 0.925 for the boundary 
layer and 0.85 for the rest of layers to 0.99 and 0.95, 
respectively. The radiative properties of cloud for this 
scheme are prescribed and are the same as those of 
CCM2. 

c. Interactive cloud radiative properties scheme: 

Fixed effective cloud droplet radius ( CWRF) 

The radiative properties of clouds are assumed to de- 
pend on their water content, the phase of the water and 
the distribution of water droplet or ice panicle sizes. In 
this scheme we attempt to compute the cloud water path 
explicitly in terms of cloud water content estimated from 
the cloud water scheme. The cloud water path, CWP, 
in units of kg m” 2 in each layer can be written as 



where p w is the cloud water density (kg m -)> C is the 
cloud fraction, q^C is the in-cloud water content, p. is 
air density, and A z is layer depth. Here, q c is the sum 
of convective cloud water and layer cloud water. For 
shortwave radiation, CWP is used to define the cloud 
extinction optical depth r f as 

r c = CWP^p + (2) 

where p and q are constants and r t is the effective radius 
of the cloud droplets (taken to be 10 p.m). The above 
formula for t c , as well as the cloud scattering and ab- 
sorption parameterizations in the following section, are 
from Slingo (1989). Here, r, will in turn affect the trans- 
missivity and reflectivity of both direct and diffuse solar 
radiation. A larger CWP for fixed r f implies an optically 
thicker cloud, and thus more reflection and less trans- 
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WRRfi CLOUD EFFECTIVE RADIUS 



CLOUD LIQUID WATER CONTENT (GM**-3) 

Fig. 1. Cloud liquid water content vs effective radius of warm 
cloud panicle (after Bower et al. 1994). 


mission of shortwave radiation. For longwave radiation, 
CWP is used to define the cloud emissivity as described 
below. 


a) — 1 - c — dr e (4) 

g = e + fr eJ (5) 

where c, d, e, and / are positive constant coefficients. 
The variation of r t will affect shortwave radiation in the 
following manner. A larger r e for fixed CWP means 
fewer droplets and a smaller effective cross section, and 
hence a smaller t c from (2). Additionally, a larger r, 
implies more absorption of radiation transmitted 
through the droplet, and hence a smaller co from (4). 
Finally, a larger r e implies more forward scattering, and 
hence a larger g from (5). 

Furthermore, in this scheme we adopt the same ab- 
sorption coefficients for long wave emissivity as those 
of Senior and Mitchell (1993), who differentiate ice 
cloud from water cloud. The cloud emissivity e, in each 
model layer k, is defined in terms of CWP (gm~ 2 ) as 

e(k) = 1 - (6) 

where the factor a is an absorption coefficient taken to 
be 0.13 m 2 g" 1 for water cloud and 0.065 m 2 g~* for 
ice cloud. Note that for the cases of CCM2, CW, and 
CWRF, a is fixed at 0.1 m 2 g'L The cloud emissivity 
is accounted for by defining an effective cloud C\ 

C f (k) = e(k)C(k) y (7) 

where C is cloud amount. Thus, larger CWP implies a 
larger emissivity, which in turn gives a larger effective 
cloud amount. 


d. Interactive cloud radiative properties scheme: 
Varying effective droplet radius for warm cloud 
only (CWRV) 


e. Interactive cloud radiative properties scheme: 
Varying effective droplet radius including ice 
cloud (CWRI) 


In addition to the dependence of cloud water path on 
cloud water content as in the case of CWRF, in this 
scheme we adopt a parameterization of the cloud droplet 
effective radius (r,) in terms of cloud liquid water con- 
tent, as suggested by Bower et al. (1994) 

r, = 1 <X)(lWC X j (3) 

where the unit of r t is pm, LWC is the liquid water 
content (in gm -3 ) predicted from the cloud water 
scheme, and N is a fixed estimate of the droplet number 
concentration, taken as 150 cm' 3 over the oceans and 
600 cm -3 over the continents (see Fig. 1). However, for 
ice cloud, the effective radius is fixed at 10 pm, which 
is the same as the prescribed value of r, for the cases 
of CCM2, CW, and CWRF. Here we are using the nu- 
merical values of Bower et al. (1994). A wide range of 
other estimates can be found in the literature. 

In addition to the dependence of the cloud extinction 
optical depth t c on r e as shown in (2), r ; is employed 
to parameterize the cloud particle single-scattering al- 
bedo o) and the asymmetry parameter g (-1 to 1) fol- 
lowing Slingo (1989); that is. 


In addition to the parameterization of warm cloud 
effective radius, in this scheme an attempt is made to 
parameterize the droplet effective radius of ice cloud in 
terms of cloud temperature from the observational study 
of Heymsfield and Platt (1984) as adapted by Suzuki et 
al. (1993) (see Fig. 2). We use 

r e = 0.717 + 61.29, (8) 

where r. is in pm and 7 is in °C. However, the upper 
and lower limits of r, are set at 40 pm and 4 pm, 
respectively, since a limitation on the size of r e is needed 
for the formulation of cloud scattering and absorption 
parameterization used in (4) and (5). In the range where 
ice and liquid coexist, the following formula (Smith 
1990) has been adopted to estimate r t . First define 



\[(T + 15)/5]- 

0 

1 - 1(775)- 


for — 15°C < T < —5°C 
for -5°C < T < 0°C. 


The resulting r e in the mixed phase can then be written 
as 
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Fig. 2. Cloud temperature vs effective radius of ice cloud particle. 
Originally from Heymsfield and Platt (1984). adapted by Suzuki et 
al. (1993). Solid line shows our parameterization. 


r t =f L r tl + ( 1-/Jr„. (9) 

where r w is from (3) and r tl is from (8) for - 15°C < 
T < 0°C. 

The five parameterizations are summarized in Table 1. 

3. Diagnostic method 

To illustrate the link between cloud cover and the 
radiation budget, the concept of cloud radiative forcing 
(Ramanathan 1987; Ramanathan et al. 1989) is used 
here. Cloud radiative forcing (CF) is defined as the ra- 
diative impact of cloud upon the top-of-atmosphere 
(TOA) radiation budget and can be expressed as 

CF = (Q - Q c ) - (F - F c ). (10) 

where Q and F represent the net incoming solar radiation 
and the outgoing longwave radiation at TOA, respec- 
tively, and the subscript C denotes clear-sky fluxes. 
Hence CF is the net cooling or heating of the earth- 
atmosphere system due to clouds. Here, CF can be sep- 
arated into its shortwave (CF SW = Q - Q c ) and long- 
wave (CF lw — F c — F) components. In general CF SW 
is negative because of the reflectance effect of cloud on 
incoming shortwave radiation; that is, Q < Q c . How- 
ever, when clouds are present, the atmospheric column 
radiates less thermal energy into space than would clear 
skies, because the effective radiating altitude is gener- 
ally greater — that is. F c > F — so that CF L%V > 0. 

To investigate the sensitivity of each scheme to 
“global warming," a ± 2 K SST perturbation is adopted 
as a surrogate climate change (Cess and Potter 1988; 
Cess et al. 1990; Zhang et al. 1994). Hereafter, 
“change” or A implies the change of a quantity arising 


Table 1. Summary of parameterizations for the five model 
versions. 



Layer cloud 
water 

Cloud optical 
thickness 

Effective radius 

CCM2 

none 

specified 

not applicable 

cw 

explicit 

specified 

not applicable 

CWRF 

explicit 

calculated 

fixed 

CWRV 

explicit 

calculated 

varying (warm cloud only) 

CWR1 

explicit 

calculated 

varying (including ice cloud) 


from this type of simulation — that is, the value of the 
quantity in the SST + 2 K case minus the value of the 
quantity in the SST — 2 K case. 

To isolate the effect of cloud feedback from that of 
other feedbacks, we adopt a method suggested by Zhang 
et al. (1994). The A F due to a climate change can be 
decomposed into the contributions from variations of 
temperature, water vapor, cloud, etc.; that is, 

A F — A F T + A F q *4- AF cl + . . . + R e , (11) 

where the subscripts T, q, and CL denote temperature, 
water vapor, and cloud, respectively, and R f is the re- 
sidual term. The shortwave components can be defined 
similarly. Therefore, the changes of cloud radiative forc- 
ings from shortwave (ACF SW ), longwave (ACF LVV ), and 
net (ACF) radiative effects, due to the simulated climate 
change, can be expressed as 


acf sw — A<2 cl 

(12) 

ACF lw = -A F cv 

(13) 

ACF = A0cl - AF C l- 

(14) 


In the computational procedure, state variables from 
the ASST = -2 K experiment are saved as (X*, k = 1, 
2, . . . ) and state variables from the ASST = F2 K 
experiment as ( X* 4- AX*, k — 1,2,...). Off-line ra- 
diation calculations are performed to obtain A F Xl 

[=F(Xj, X., . . . , X* + AX*, . . . ) - F(X„ X. X*, 

. . . )] and AG*, [ = Q(Xj, X : X* + AX*, . . . ) - 

g(X„ X : , . . . , X*, . . . )] at each grid point. A total of 
240 fields, obtained by sampling instantaneously every 
6 simulated hours over a total of 60 simulated days, are 
used in this study, in view of the calculated magnitude 
of differences in radiative fluxes between the simulated 
atmospheric states sampled hourly and the instantaneous 
fields. The partial feedbacks from cloud are attributed 
to variations of both cloud fraction and cloud micro- 
physical properties, such as cloud water path and cloud 
droplet effective radius. Our sampling strategy is meth- 
od 4 of Zhang et al. (1994). 

This diagnostic method of off-line radiation calcu- 
lation has also been applied to the investigation of the 
partial cloud radiative forcing of high, middle, and low 
cloud of the simulated present climate, which will be 
discussed below. 




2484 


JOURNAL OF CLIMATE 


VOLUME 10 


Table 2. Radiation budget and cloud radiative forcing (W m~ 2 ) 
for ERBE observations and the five model versions. See text for 
definitions. 



ERBE 

CCM2 

CW 

CWRF 

CWRV 

CWRI 

Q 

231.2 

243.3 

236.6 

232.5 

222.8 

235.6 

Qc 

277.4 

288.1 

287.9 

287.9 

287.7 

287.8 

F 

238.3 

246.1 

233.0 

240.3 

243.2 

244.6 

F c 

266.6 

275.3 

273.4 

273.7 

272.6 

273.4 

CF sw 

-46.2 

-44.8 

-51.3 

-55.3 

-64.9 

-52.3 

CF lw 

28.3 

29.2 

40.4 

33.4 

29.4 

28.8 

CF 

-17.9 

-15.6 

-10.9 

-21.9 

-35.5 

-23.4 


4. Model simulations 

a. Present climate 

In Table 2, the global mean simulated and observed 
cloud radiative forcings are summarized. All the cases 
show negative net cloud forcing; that is, cloud has a 
cooling effect on the earth-atmosphere system. The dis- 
crepancies between ERBE observations and the model 
results arise from both cloudy-sky and clear-sky fluxes. 
For the resulting global mean cloud radiative forcing, 
the reasonable agreement between ERBE and CCM2 is 
noteworthy. 

The simulated radiation budget may be explained in 
terms of the variables listed in Table 3, which sum- 
marizes the global mean values of three model- 
generated variables: total fractional cloud coverage, ver- 
tically integrated cloud water path, and effective cloud 
droplet radius. 

For the estimation of the global mean effective radius, 
only grid points where the fraction is not zero are count- 
ed. For the cases without interactive cloud radiative 
properties (CCM2 and CW), the effects of cloud on the 
radiation budget and the consequent cloud radiative 
forcing are mainly attributable to the cloud fractional 
coverage. If we compare CCM2 and CW, the smaller 
cloud fraction of CCM2, compared to that of CW, sim- 
ply explains the larger values of both Q and F of CCM2, 
compared to those of CW. However, for the cases with 
variable cloud radiative properties (CWRF, CWRV, and 
CWRI), we must invoke additional variables to interpret 
the results. The differences in Q for the cases of CWRF, 
CWRV, and CWRI. where the total cloud fractions are 
not much different, may be explained by the differences 
of r € through (2), (4), and (5). For the case of CWRV, 
the smallest r t results in the smallest Q> and for the case 
of CWRI, the opposite is true; that is, a larger r e implies 
a smaller optical depth from (2), and a smaller particle 
scattering albedo from (4), which leads to larger short- 
wave fluxes. For longwave radiation, we notice that F 
depends on cloud water path, if cloud fractional cov- 
erages are comparable; that is, a larger cloud water path 
implies a larger cloud emissivity from (6), and conse- 
quently a smaller F. 

Figure 3 shows the zonal mean distribution of total 
cloud fraction for each case. CCM2 produces the small- 


Table 3. Global mean total cloud, vertically integrated cloud water 
path, and effective cloud droplet radius for the five model versions. 



CCM2 

cw 

CWRF 

CWRV 

CWRI 

Cloud (%) 

55.4 

63.8 

66.1 

66.1 

65.5 

Cloud water path (gm' 2 ) 

403.5 

434.4 

193.6 

189.1 

185.8 

Effective radius (ptm) 

10.00 

10.00 

10.00 

9.67 

11.31 


est cloud fraction at virtually every latitude. Figure 4 
shows the zonal mean distribution of shortwave and 
longwave radiation and the associated cloud radiative 
forcings for the simulated fields and the ERBE data. 
The zonal mean distribution of the resulting net radiative 
forcing CF is shown in Fig. 5. It is obvious from these 
figures that the reasonable agreement in global mean 
cloud radiative forcing (Table 2) between ERBE and 
CCM2 is largely due to canceling effects at different 
latitudes. For the cases of CWRF, CWRV, and CWRI, 
in which cloud radiative properties are explicitly param- 
eterized, a sizeable improvement over CCM2 and CW, 
in which cloud radiative properties are specified, is seen 
in Q and F in the middle latitudes of the Northern Hemi- 
sphere. 

Figures 6, 7, and 8 depict the height-latitude cross 
sections of cloud fraction in percentages, cloud water 
path in gm -2 , and cloud droplet effective radius in ft m, 
respectively. In Fig. 6 we show only the cases of CCM2 
and CWRF since differences in the distributions of 
cloud fraction among CW, CWRF, CWRV, and CWRI 
are barely distinguishable. This is true for Fig. 7 as well, 
where we show only the cases of CCM2 and CWRI. 
The cloud water paths of CCM2 and CW, in which those 
are specified in terms of height and layer thickness, are 
virtually identical to one another, and those of CWRF, 
CWRV, and CWRI, in which the same formula for cloud 
water path is applied, are also almost indistinguishable 


TOTRL CLOUD FRRCT I ON 



Fig. 3. Zonally averaged total cloud fractions in percents from 
60-day perpetual July simulations. 
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Q 4 QC Q-QC 



F & FC FC-F 



ERBE 

COM2 

cw 

CWRF 

CWRV 

EwrT 


Fig. 4. Zonally averaged shortwave (0 and longwave (F) radiation fluxes at the top of the 
atmosphere, their clear-sky components ( Q c and F c ), and the resulting shortwave and longwave 
cloud radiative forcings ( Q — Q c and F c — F, respectively) from 60-day perpetual July simu- 
lations and the ERBE observational data for July 1985-88. 


from one another. The distributions of effective radius 
of CWRV and CWRI are shown in Fig. 8. bearing in 
mind that CCM2, CW, and CWRF have a fixed effective 
radius of 10 /zm. In Fig. 8, we easily notice the effect 
of the ice cloud effective radius parameterization of 
CWRI; in the regions of warm cloud, there is little dif- 
ference between the two models in the effective radius, 
while for the regions of ice cloud, the effective radius 
of CWRI is bigger than that of CWRV due to the 
temperature-dependent ice cloud effective radius param- 
eterization of (8). 

To investigate the height-dependence of cloud radi- 
ative forcing, we utilize the diagnostic method of off- 
line radiation calculation. For this purpose, we define 
clouds above 400 mb as high cloud, between 400 and 
700 mb as middle cloud, and below 700 mb as low 
cloud. Figure 9 summarizes the results. For the short- 
wave forcing, the total forcing is mostly from the op- 
tically thick (large CWP) low cloud, while for the long- 
wave forcing, it is mostly from the high cloud, as noted 
in other studies (e.g., Senior and Mitchell 1993). In 
terms of net cloud radiative forcing, the effect of high 
cloud is positive, and that of middle and low clouds is 


negative, giving a net negative cloud forcing. Note that 
there is a slight difference in the estimations of total 
cloud forcing between Fig. 9 and Table 2, since the 
former is estimated from the off-line radiation calcu- 
lation, while that of Table 2 is from the model-simulated 
fields. Also note that the total forcing is not the same 
as the sum of the forcings from high, middle, and low 
cloud, since the assumption of random overlapping in 
the vertical direction emphasizes the highly nonlinear 
relationship between the total radiative forcing and the 
separate forcings, as also pointed out by Zhang et al. 
(1994). 

b. Sea surface temperature fluctuation experiments 

Table 4 summarizes the global means of feedbacks 
including partial contributions from cloud, water vapor, 
and temperature for each case. “FULL” denotes the 
estimation from the hourly averaged simulation, while 
“SUM” denotes the sum of partial terms (see also Fig. 
10). Figure 11 shows the zonal mean changes in radi- 
ative fluxes due to the ASST climate change. The dif- 
ferences between the sum of the partial terms and the 
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NET CLOUD RADIATIVE FORCING 



FlC. 5. Zonally averaged net cloud radiative forcings from 60-day 
perpetual July simulations and the ERBE observational data for July 
1985-88. 


simulated variation are indeed small, indicating that the 
variations of radiative fluxes are almost linear for this 
climate change of about AT = 4 K. For the changes in 
shortwave forcing, most of the variation is due to the 
changes in cloud, but the contribution from changes in 
water vapor, resulting from the increase of absolute hu- 
midity in the warmer climate, are also noticeable. There 
is no contribution in A Q from the changes in temper- 
ature, as expected. For the changes in longwave forcing, 
the water vapor greenhouse effect is apparent in terms 
of the reduction in F q . The sum of AF r and A F qy which 
expresses the net effect of temperature and water vapor 
on the longwave flux, is positive; that is, there is a net 
radiative loss of energy due to feedbacks involving these 
two factors. The thin solid lines in Fig. 1 1 represent the 
zonal mean distribution of cloud radiative forcing. Table 
5 denotes the global mean cloud radiative forcing based 
on data from Table 4, and for comparison purposes, the 
estimates from the method using (10) are denoted by 
(old). It is noteworthy that the signs of cloud radiative 
feedbacks are different between the schemes with spec- 
ified cloud radiative properties — that is, CCM2 and 
CW — and the schemes with interactive radiative prop- 
erties — that is, CWRF, CWRV, and CWRI. For the cases 
of CCM2, CW, and CWRI, the net cloud radiative feed- 
back is negative — that is, cloud moderates the global 
warming, while for the cases of CWRF and CWRV, the 
opposite is true. 

A useful measure of the cloud radiative feedback is 
the quantity ACRF/G, which is the change in cloud 
radiative forcing normalized by the direct radiative forc- 
ing (e.g., Cess et al. 1996). When ACRF/G is unity, the 
direct radiative forcing is amplified by a factor of 2, 



CCM2 


CLIRF 


Fig. 6. Height-latitude cross section of cloud amount for CCM2 
and CWRF from 60-day perpetual July simulations. Contours every 
5%; regions of over 25% are stippled. 


that is, a twofold positive feedback. In terms of the 
variables in Table 4, ACRF/G is equivalent to -A (Q - 
F) C JA(Q — F ) FULL . This quantity is shown in Table 4, 
where it can be seen to vary from —0.18 to 0.32. Thus, 
the range of cloud radiation parameterizations explored 
in these experiments yields feedbacks ranging from 
slightly negative to slightly positive on a global basis. 
Of course, the local effects of these feedbacks may be 
more significant than the global mean values, which are 
the focus of the present study. 

Table 6 is similar to Table 3 except for the changes 
in the ASST experiments. The positive ACF SW for the 
cases of CCM2 and CW is simply related to the reduc- 
tion of total cloud fraction. Similarly, the negative 
ACF lw for the cases of CCM2 and CW is also related 
to the reduction of total effective cloud fraction. How- 
ever, for the cases of CWRF, CWRV, and CWRI, the 
increase in CWP results in a positive ACF LW even 
though the total cloud fraction decreases. For shortwave 
radiation, the negative ACF SW for the case of CWRF is 
due to the increase in CWP, which dominates over the 
effect of a decrease in total cloud fraction, while the 
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latitude latitude 

Fig. 7. Height-latitude cross section of cloud water path for CCM2 Fig. 8. Height-latitude cross section of cloud droplet ettective 
and CWRI from 60-day perpetual July simulations. Contours every radius for CWRV and CWRI from 60-day perpetual July simulations. 

10 girr 2 ; regions of over 50 gm~ 2 are stippled. Contours every 1 fim; regions of over 10 /am are stippled. 


negative ACF SW for the cases of CWRV and CWRI is 
due to the combined effects of an increase in CWP and 
a decrease in r e . 

The literature contains conflicting estimates of the 
temperature dependence of cloud optical thickness. For 
a recent discussion, see Del Genio et al. (1996). 

Figures 12, 13, and 14 are similar to Figs. 6, 7, and 
8 except for the changes in the A SST experiments. In 
Fig. 12, the increase in high cloud at low latitudes and 
the decrease in middle-level cloud are common to all 
of the models, even though here we present only the 
cases of CCM2 and CWRF. The maximum in high cloud 
at low latitudes moves upward as the height of the tro- 
popause increases, as found in previous numerical stud- 
ies (e.g., Mitchell and Ingram 1992). In Fig. 13, the 
change in cloud water path of CCM2 is indeed negli- 
gible compared with that of CWRI, as also noted in 
Table 6. In Fig. 14, there is a large reduction in the 
cloud effective radius for the case of CWRI at middle 
latitudes, which is due to a decrease in the temperature- 
dependent effective radius of ice and mixed cloud — that 


is, to a change of state from ice (larger r e ) to liquid 
(smaller r,). 

Figure 15 is similar to Fig. 9. but again depicts the 
changes. The decrease in longwave forcing for the cases 
of CCM2 and CW is mainly attributable to a decrease 
in high cloud fraction. However, for the shortwave feed- 
back of CCM2, it is noteworthy that an increase in low 
cloud water path, even though its magnitude is small, 
results in a negative low cloud shortwave feedback. The 
increase in high cloud fraction for the cases of CWRF, 
CWRV, and CWRI results in a positive high cloud net 
radiative feedback. The smaller negative high cloud 
feedback in CWRI compared to CWRV may be ascribed 
to the combined effects of a smaller increase in high 
cloud fraction and a larger increase in cloud droplet 
effective radius. However, the decrease in middle cloud 
effective radius does not have much impact on the re- 
sulting middle cloud shortwave forcing. The shortwave 
cloud radiative feedback due to effective radius is much 
smaller than that of cloud water path. Even though the 
changes in low cloud amount and low cloud water path 
are apparent, the resulting change in low cloud long- 
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Fig. 9. Globally averaged total cloud fraction, vertically integrated cloud water path, cloud 
droplet effective radius, shortwave, longwave, and net cloud radiative forcings (denoted by O). 
Also shown are results from separate high, middle, and low cloud, denoted by A, +, and X, 
respectively. All results are from 60-day perpetual July simulations. 


Table 4. Diagnostic analysis of feedbacks (W m _: ). See text for 
definitions. 




CCM2 

CW 

CWRF 

CWRV 

CWRI 



3.60 

5.01 

-1.15 

-2.03 

-3.04 

A£?sum 


3.24 

4.70 

-1.53 

-2.42 

-3.44 

A<2r 


0.00 

0.00 

0.00 

0.00 

0.00 

as. 


0.83 

0.68 

0.68 

0.83 

0.71 

A0CL 


2.41 

4.02 

-2.20 

-3.25 

-4.15 

A/W 


14.82 

13.44 

4.95 

5.69 

5.83 

AF SUM 


15.21 

13.98 

4.67 

5.34 

5.51 

AF r 


18.19 

15.88 

17.30 

17.12 

17.11 

AF, 


-7.39 

-6.06 

-8.46 

-8.32 

-8.29 

af cl 


4.40 

4.16 

-4.17 

-3.45 

-3.31 

A(G - 


-11.22 

-8.43 

-6.10 

-7.72 

-8.87 

A«2 - 

f) S LM 

-11.97 

-9.28 

-6.20 

-7.76 

-8.96 

A (G - 

F) r 

-18.19 

-15.88 

-17.30 

-17.12 

-17.11 

A«2 - 

n. 

8.22 

6.74 

9.13 

9.16 

8.99 

A(Q - 

F ) CL 

-1.99 

-0.14 

1.96 

0.20 

-0.84 

ACRF IG 

-0.18 

-0.02 

0.32 

0.03 

-0.09 


wave forcing is negligible. For the net change in cloud 
radiative forcings, we notice that high cloud produces 
a positive feedback, while low and middle clouds pro- 
duce a negative feedback. 

5. Observational validation 

The preceding results demonstrate clearly that GCM 
simulations of climate sensitivity show a strong depen- 
dence on the way in which cloud radiation processes 
are parameterized. However, these model results cannot 
settle the question of which parameterizations are most 
realistic. Only a thorough comparison with observations 
can definitively show that one algorithm is superior to 
another. For example, it is not correct that the most 
physically comprehensive treatment is necessarily the 
best. It is one thing to design a scheme in which, for 
instance, radiative properties depend on a computed 
cloud water budget, and it is quite another to do so 
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Fig. 10. Changes in globally averaged radiation (O) and its partition into variations of 
cloud (A), water vapor ( + ), and temperature (X) due to 4SST = ±2 K experiments con- 
sisting of 60-day perpetual July simulations. 


realistically. Simply to assert that a given physical pro- 
cess is “included" or “taken into account" conveys 
nothing about whether the result is actually an improve- 
ment on a less ambitious scheme. For example, a simple 
prescription of cloud radiative properties may well be 
more appropriate than an overly complicated attempt to 
compute them based on questionable assumptions about 
cloud microphysics. 

How then should one go about assessing the merits 
of parameterizations such as those described in section 
2? There are many possible routes, and a comprehensive 
assessment is still a goal for the future. Nevertheless, 
recent progress has made it possible to begin the long- 
needed task of testing GCM algorithms empirically. In 
the present study, we have chosen to illustrate the po- 
tential of recent observational and theoretical advances 
by using a single-column model (SCM) in conjunction 
with measurements from two field programs, ARM and 
TOGA COARE. In this section, we describe the SCM 
and the field programs briefly, and we show results from 
using the SCM to validate the parameterizations used 
in our GCM tests against observations from the two 
field programs. Both these field programs provide in- 
tensive measurements of an atmospheric volume com- 
parable to that overlying a single GCM horizontal grid 


cell. One program (ARM) provides data from a mid- 
continent, midlatitude region (the central plains of the 
United States). The other program (TOGA COARE) 
took place in the western tropical Pacific Ocean. Thus, 
we have independent data from two sites, w'hich are 
dramatically different meteorologically. The theoretical 
tool for using these data to test GCM parameterizations 
is the SCM. 

a . SCM description 

The SCM is a computationally efficient and econom- 
ical one-dimensional (vertical) model, resembling a sin- 
gle column of a general circulation model (GCM). The 
SCM is applied at a specific site having a horizontal 
extent typical of a GCM grid cell. Since the model is 
one-dimensional, the advective terms in the budget 
equations must be specified from observations or op- 
erational numerical weather prediction analyses. 

The single-column mode! is diagnostic rather than 
prognostic. Its input is an initial state, plus the time- 
dependent advection terms in the conservation equa- 
tions, provided at all model layers. Its output is a com- 
plete heat and water budget for the study site, including 
temperature and moisture profiles, clouds and their ra- 
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Fig. 11. Changes in zonally averaged variation of radiation and its partition into variations 
of temperature, water vapor, and cloud due to ASST = ±2 K experiments consisting of 60-day 
perpetual July simulations; “FULL” (thick solid line), "SUM” (thick dotted line), "cloud’ 
(thin solid lines), “water vapor” (thin dotted line), and "temperature” (thin dashed line). 


diative properties, diabatic heating terms, surface energy 
balance components, and hydrologic cycle elements, all 
specified as functions of time. For a more complete 
discussion of SCMs, the reader is referred to Randall 
et al. (1996). 

The present configuration of the single-column model 
has 16 vertical layers in the atmosphere and uses a time 
step of 7.5 min. A diumally varying solar signal de- 
pendent on the latitude and time of year is applied at 
the top of the model atmosphere using a solar constant 
of 1370 W nr 2 . 

The SCM can be run in a “fully interactive” mode, 


Table 5. Cloud radiative feedbacks (W m : ) for the five model 
versions. See text for definitions. 



CCM2 

CW 

CWRF 

CWRV 

CWRI 

acf sw 

2.41 

4.02 

-2.20 

-3.25 

-4.15 

acf lw 

-4.40 

-4.16 

4.17 

3.45 

3.31 

ACF 

-1.99 

-0.14 

1.96 

0.20 

-0.84 

ACF SW (old) 

2.97 

4.57 

-1.75 

-2.48 

-3.50 

ACF LW (old) 

-6.08 

-5.36 

2.45 

1.86 

1.54 

ACF (old) 

-3.11 

-0.79 

0.70 

-0.62 

-1.96 


in which after being initialized, the SCM determines the 
temperature and humidity profiles at each time step us- 
ing the heating and moistening rates calculated by the 
various model parameterizations. The main disadvan- 
tage of this mode is that accumulated model errors in 
the temperature and humidity profiles may affect the 
performance of the model parameterizations. Alterna- 
tively. the SCM can be run in the “semi-implicit” mode 
where the temperature and humidity profiles are spec- 
ified from observational data at each time step. A draw- 
back of this technique is that it may distort the atmo- 
spheric profiles that could affect significantly the heating 

Table 6. Changes in global mean total cloud, vertically integrated 
cloud water path, and effective cloud droplet radius for the five model 
versions. See text for definitions. 

CC M2 CW CWRF CWRV CWRI 

Ac loud ( a c) -1.33 -1.59 -0.50 -1.08 -1.04 

Acloud water 

path (gm" 2 ) 0.000061 0.000012 66.33 61.11 66.46 

^effective 

radius (^m) 0.000 0.000 0.000 -0.032 -0.143 
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LRTITUCE 

Fig. 12. Height-latitude cross section of equilibrium changes in 
cloud amount for CCM2 and CWRF from 60-day perpetual July ±2 
K SST simulations. Contours every 1 9c\ regions of increase are stip- 
pled. 


LATITUDE 

Fig. 13. Height-latitude cross section of equilibrium changes in 
cloud water path of CCM2 and CWRI from 60-day perpetual July 
±2 K SST simulations. Contours every 0.1 gm 2 and 10 grrr 2 of 
CCM2 and CWRI, respectively; regions of reduction are stippled. 


and moistening rates produced by the cumulus convec- 
tion parameterization. 

A compromise between these two modes is to operate 
the SCM in “relaxation mode’’ in which the temperature 
and humidity profiles are relaxed to the observed pro- 
files using a specified time constant, (3. A value of /3 = 
oo represents the fully interactive mode, while a value 
of j8 = A/, where A t is the model time step (7.5 min 
in this study), represents the semi-implicit mode. In the 
SCM experiments described in this paper, a time con- 
stant of P = 24 h is used. The choice of 24 h helps to 
ensure that model errors in the temperature and humidity 
profiles remain small without seriously distorting the 
atmospheric profiles. 

The SCM contains a full set of parameterizations of 
subgrid physical processes that are normally found in a 
modem GCM. In the SCM experiments that follow, the 
parameterizations used to determine the cloud properties 
are the same as in the five schemes CCM2, CW, CWRF, 
CWRV, and CWRI. The SCM does employ some phys- 
ical parameterizations that differ from those used in 


CCM2. These include the cumulus convection param- 
eterization of Zhang and McFarlane (1995), the solar 
radiation parameterization of Fouquart and Bonnel 
(1980), and the longwave radiation parameterization of 
Morcrette (1990). 

SCM integrations were performed at two different 
locations. The first group used data collected during 
TOGA COARE located in the western tropical Pacific. 
The second group of integrations were located over the 
U.S. southern Great Plains and used data collected from 
the ARM program. Both sets of SCM integrations use 
in situ data to calculate the necessary horizontal ad- 
vective forcing terms. The main advantage in specifying 
the forcing terms using in situ data as opposed to three- 
dimensional model analysis products is that the latter may 
contain biases from the physical parameterizations con- 
tained in the three-dimensional model. Each group of SCM 
integrations consisted of five model runs with the same 
cloud parameterizations as in the CCM2, CW, CWRF, 
CWRV, and CWRI schemes described in section 2. 

Satellite measurements of outgoing longwave radia- 
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humidity, and horizontal winds were produced every 6 
h at each of these four sites during TOGA COARE 
throughout the 4-month IOP period of 1 November 
1992-28 February 1993. The original radiosonde mea- 
surements were quality controlled (for details, see Miller 
and Riddle 1994) and then interpolated onto evenly 
spaced vertical levels with constant pressure increments 
of 5 mb. Surface observations of temperature, humidity, 
wind, precipitation, and radiative fluxes were also ob- 
tained at each site and stored as 1-min averages. 

Sounding data from these four sites were used to 
produce slab averages of temperature, humidity, and 
horizontal winds for the IFA column. The horizontal 
divergence at each pressure level was calculated by eval- 
uating the line integral around the IFA perimeter: 


1 

V V = - 
A 


V*n dl. 


where V is the horizontal velocity, n is the normal vector 
relative to the IFA perimeter, and A is the area enclosed 
by the IFA perimeter. Vertical velocities were calculated 
from the horizontal divergences using the method of 
O’Brien (1970). 

The horizontal advective tendency of temperature was 
calculated using 


v-vr 



rv-n dl 


T 


Vn dll 


Fig. 14. Height-latitude cross section of equilibrium changes in 
cloud droplet effective radius for CWRV and CWRI from 60-day 
perpetual July ±2 K SST simulations. Contours every 0.5 /i,m; 
regions of reduction are stippled. 


tion and cloud fraction together with surface measure- 
ments of downwelling shortwave radiation are used to 
validate the model results from each of the five cloud 
parameterization schemes. By using these three vari- 
ables for comparison, changes in overall cloud fraction, 
cloud optical thickness, and cloud-top height are con- 
sidered. The in situ data used to calculate the forcing 
terms and the data used to validate the SCM results are 
discussed below for each location. 

b. SCM experiments — Western tropical pacific 

The timing of these SCM integrations coincides with 
the 4-month Intensive Operations Phase (IOP) of TOGA 
COARE. The SCM was applied in the atmospheric col- 
umn directly over the TOGA COARE Intensive Flux 
Array (IFA). Two land-based sites and two research ves- 
sels deployed at stationary sites defined the perimeter 
of the TOGA COARE IFA, which covers an area of 
approximately 240000 km : . 

Atmospheric soundings consisting of temperature, 


where T is the temperature and the overbar denotes a 
spatial average. A similar equation is used for the hu- 
midity tendency. 

Due to missing and errant data there were several 
significant gaps in the sounding measurements. As a 
result, the SCM experiments were limited to three time 
periods for which near-continuous forcing data was 
available. These three periods (10 November 1992-11 
December 1992, 18 December 1992-22 January 1993, 
and 31 January 1993-19 February 1993) contain a total 
of 88 days. 

Outgoing longwave radiative flux and fractional cloud 
amount were obtained from the International Satellite 
Cloud Climatology Project (ISCCP) C 1 flux dataset cre- 
ated in support of TOGA COARE. The dataset is pre- 
sented with a temporal resolution of 3 h and a spatial 
resolution of 2.5° X 2.5°. Downwelling surface short- 
wave radiative fluxes were obtained from a total of sev- 
en surface sites within the IFA. These include the four 
ISS sites defining the perimeter of the IFA, two ships 
(R/V Moana Wave and R/V PRC#5), and the Woods 
Hole Oceanographic Institution (WHOI) IMET buoy. 
These seven surface measurements were averaged to 
produce a single time series of downwelling surface 
shortwave radiation representative of the entire IFA re- 
gion. 
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Fig. 15. As in Fig. 9 but for the changes due to JSST = ±2 K experiments. 


c. SCM experiments — U.S. southern Great Plains 

This set of SCM integrations utilizes data collected 
at the ARM southern Great Plains (SGP) site during the 
21-day fall 1994 IOP (25 October 1994-14 November 
1994). During this IOP, radiosondes were launched si- 
multaneously from five different locations every 3 h. 
These five locations consist of four boundary sites and 
one central site. The size of the array defined by these 
boundary sites is approximately 250 km X 250 km. The 
radiosondes measure the u- and i>wind components, 
temperature, and humidity as a function of pressure. 
Collocated with the radiosonde launch sites are National 
Weather Service 404-MHz radar wind profilers. 

Temperature and humidity data from each radiosonde 
were averaged over 10-mb increments extending from 
the surface to 100 mb to form “layer" averages. Ob- 
jective analysis techniques were used to combine the 
radiosonde wind data and the wind profiler data into 
layer averages. The wind divergence and horizontal ad- 
vective tendencies were then calculated by differencing 
the appropriate layer averages of the four boundary 
sites. Linear interpolation was used (both spatially and 
temporally) to fill in any missing data. 

Outgoing longwave radiation and fractional cloud 


amount were derived from GOES-7 measurements 
(Minnis et al. 1995) and obtained from the NASA/Lang- 
ley Research Center. A limitation of this dataset is that, 
due to a lack of cloud optical depth information during 
nighttime scenes, the nighttime cloud fraction should 
be regarded with caution. Daytime values of cloud frac- 
tion should not be affected. In the data comparisons that 
follow the total cloud fraction is used in the plots. How- 
ever, the daytime cloud fractions from both the SCM 
and the GOES-7 products are noted in the text. 

Downwelling shortwave radiation measurements 
were obtained from a total of 31 surface sites within 
the ARM SGP site boundary. These included measure- 
ments from the ARM central site facility, 2 sites from 
the Kansas State University Mesonetwork, and 28 sites 
from the Oklahoma Mesonetwork (Brock et al. 1995). 
The data from these 31 sites were averaged to produce 
a single time series of downwelling surface shortwave 
radiation representative of the ARM SGP site. 

d. SCM results 

The model-produced outgoing longwave radiation, 
downwelling surface shortwave radiation, and cloud 
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Fig. 16. Comparison between model and observed mean outgoing 
longwave radiation for each of the five cloud schemes. The results 
from the SCM experiments over the western tropical Pacific are 
shown in the top panel while the bottom panel shows the results from 
the U.S. southern Great Plains. 

fraction from each SCM integration (a total of 10) were 
averaged over the length of the model run and then 
compared to the corresponding observed means. Figure 
16 shows the comparison between model and observed 
outgoing longwave radiation for each of the five cloud 
schemes. At both the TOGA COARE site and the ARM 
SGP site the incorporation of prognostic cloud liquid 
water significantly improves the model outgoing long- 
wave radiation results (CCM2 vs CW). When the cloud 
radiative properties are computed as a function of cloud 
liquid water and/or cloud droplet radius (CWRF, CWRV, 
or CWRI), further improvement is obtained. However, 
the improvement is not as great as that between CCM2 
and CW. 

The modeled and observed mean cloud fractions are 
shown in Fig. 17. Again, the model results at both sites 
show a significant improvement with the incorporation 
of prognostic cloud liquid water. However, very little 
difference in the modeled mean cloud fraction is seen 
when cloud radiative properties are computed as a func- 
tion of cloud liquid water and/or cloud droplet radius 
(CW vs CWRF, CWRV, or CWRI). The observed cloud 
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Fig. 17. As in Fig. 16 but for mean cloud fraction. 


fraction at the ARM SGP site shown in Fig. 17 is 0.55 
and was calculated using both day and nighttime values. 
If only daytime values are considered, the observed 
cloud fraction increases to 0.61. The modeled daytime 
mean cloud fraction varied less than 5% from the values 
shown in Fig. 17. Thus, the limitations of the GOES-7 
cloud products discussed above do not appear to alter 
the results of these comparisons. 

Figure 18 shows the modeled and observed mean 
downwelling surface shortwave radiation. As with the 
outgoing longwave radiation and the cloud fraction, a 
notable improvement in the modeled downwelling sur- 
face shortwave radiation occurs with the inclusion of 
prognostic cloud liquid water. When cloud radiative 
properties are determined from cloud liquid water and/or 
cloud droplet radius, an additional, albeit smaller, im- 
provement is obtained. 

Thus, the inclusion of prognostic cloud water does 
seem to produce a notable improvement in the realism 
of the parameterizations. as measured by SCM diagnosis 
of observations from TOGA COARE and ARM. Of 
course, the observations themselves have their own 
shortcomings, as we have noted. Nevertheless, it seems 
fair to conclude that empirical validation of the tested 
algorithms demonstrates the value of including cloud 
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Fig. 18. As in Fig. 16 but for mean downwelling surface short- 
wave radiation. 


water as a prognostic variable. Furthermore, the obser- 
vational evidence suggests that deriving cloud radiative 
properties from cloud water content and microphysical 
characteristics is a promising route to further improve- 
ment. 

6. Conclusions 

The cloud cover and radiation fields for simulations 
of the present climate are strongly sensitive to the cloud 
radiation scheme adopted. Companson with satellite ob- 
servations shows that improvements in all the algo- 
rithms are still needed. In the warming experiments, we 
find either a negative or a positive net cloud feedback 
depending on which schemes we chose. However, for 
all the cases with interactive cloud radiative schemes, 
we find a negative shortwave feedback (i.e., cloud mod- 
erates the warming through its effect on shortwave ra- 
diation) and a positive longwave feedback (i.e., cloud 
amplifies the warming through its effect on longwave 
radiation), while the opposite is true for the schemes 
with specified cloud radiative properties. The increase 
in cloud water content in the warmer climate leads to 
optically thicker middle and low clouds and in turn to 
negative shortwave feedbacks for the interactive cloud 


radiative schemes, while the decrease in cloud amount 
simply produces a positive shortwave feedback for the 
other schemes. For the longwave feedbacks, the de- 
crease in high effective cloudiness for the schemes with- 
out interactive radiative properties leads to a negative 
feedback, while the opposite is true for the schemes with 
interactive cloud radiative properties. 

Our results reinforce a major conclusion of Senior 
and Mitchell (1993), who showed that variations in the 
cloud radiation algorithm of their GCM could account 
for much of the typical factor-of-three difference be- 
tween GCMs (in global sensitivity of surface temper- 
ature to greenhouse gas increases) that has been revealed 
by extensive model intercomparisons. Like they and 
others, we find significant differences in the magnitudes 
and even the signs of cloud radiation feedbacks, de- 
pending on our choice of parameterizations. It is some- 
what unsettling that the results of a modem atmospheric 
GCM, which may contain several tens of thousands of 
lines of code, can be affected so much by making what 
seem to be reasonable changes in only a few of those 
lines. This sensitivity highlights the need for a thorough 
understanding of the fundamental physics of clouds and 
their effects on radiation. 

Unfortunately, it is not known which cloud radiation 
parameterizations are most realistic, or even if any cur- 
rent scheme captures the essential feedback processes 
of actual clouds. Furthermore, other GCM groups have 
obtained different results by trying other ways of in- 
corporating cloud microphysical processes and their ra- 
diative interaction. There is thus a clear need to inter- 
compare these approaches with one another, and. even 
more importantly, to validate them against observations, 
so as to evaluate the strengths and weaknesses of each. 
For too long, we have had too many theoretically plau- 
sible parameterizations and too little research on testing 
them empirically. By using the single-column model in 
conjunction with observations from ARM and TOGA 
COARE, we have carried out a partial validation of 
several cloud radiation algorithms. The inclusion of 
prognostic cloud water produces a notable improvement 
in the realism of the parameterizations. as measured by 
observations. We conclude that empirical validation of 
the tested algorithms demonstrates the value of includ- 
ing cloud water as a prognostic variable. Furthermore, 
the observational evidence suggests that deriving cloud 
radiative properties from cloud water content and mi- 
crophysical characteristics is a promising route to fur- 
ther improvement. 
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